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ASYMPTOTIC  DISTRIBUTION  OF  THE  SHAPIRO-WILK.  W  FOR  TESTING  FOR  NORMALITY 

By 

J.  R.  Leslie,  M.  A.  Stephens  and  S.  Fotopoulos 

1.  Introduction. 

A  popular  test  for  the  normality  of  a  random  sample  is  based  on  the 

Shapiro-Wilk  statistic  W.  This  statistic,  which  was  presented  in  Shapiro 

and  Wilk  (1965),  is  the  ratio  of  the  square  of  the  BLUE  of  a  to  the 

2 

sample  variance,  where  a  is  the  variance  of  the  normal  population  from 

which  the  samule  is  assumed,  under  the  null  hypothesis,  to  have  been 

1/2 

drawn.  For  convenience  we  shall  work  with  W  which  has  the  form 

w1/2  =  X ’ V^m/ { 7  (X.-X ) 2m* V^V" 1 ' 2  , 

where  X  =  (X, , . . . ,X  ) ' ,  X,  <  X0  <  . . .  <  X  ,  is  the  vector  of  order 
—  Ini/  n 

statistics  from  the  sample,  X  is  the  sample  mean,  and  m  is  the  mean 

vector  and  Vq  the  covariance  matrix  of  standard  normal  order  statistics. 

1/2 

As  W  is  location  and  scale  invariant  we  can  assume  from  henceforth 

that  X. , . . . ,X  ,  are  order  statistics  for  a  sample  from  a  N(0,1) 
l  n 

population . 

A  number  of  authors  (for  example,  Sarkadi  (1975),  (1977)  and  Gregory 

(1977))  have  (correctly)  guessed  at  the  form  of  the  asymptotic  distribution 

for  W  as  well  as  predicting  that  the  test  should  be  consistent.  However 

no  rigorous  proofs  have  been  possible  due  to  the  presence  of  V^. 

Neither  V'  nor  can  be  found  explicitly  and  until  recently  no 

0  0 


1 


reasonably  accurate  asymptotic  approximation  for  Vq  was  available.  A 
paper  by  one  of  the  authors  (Leslie  1984)  has  now  remedied  the  situation; 
in  that  paper  can  be  found  an  approximation  for  Vq  together  with  a  number 
of  asymptotic  properties  of  Vq,  one  of  which  is  of  particular  importance 
to  this  work.  It  states  that  m  is  approximately  an  eigenvector  of  V ^ 
in  the  following  sense: 

(1)  //  Vq1!- 2m//  <  C(logn)_1/2  , 

2  2 

where  C  is  a  constant  independent  of  n,  and  //  b //  =  Ib^  for 

b  =  (b^,...,b  )'.  This  latter  result  formalises  a  similar  one  appearing 
in  Stephens  (1975). 

The  asymptotic  distribution  of  W,  after  appropriate  normalizing, 
has  been  assumed  to  be  the  same  as  that  of  the  De  Wet  and  Venter  (1972) 
statistic 

W*  =  r2(X,H)  ; 

here  r(X,Y)  is  the  sample  correlation  coefficient  between  X  and  Y, 

th 

H  is  the  n*  1  vector  whose  i  element  is  $  (i/(n+l))  and  $  (•) 

is  the  inverse  function  for  the  standard  normal  distirbution  function 
$(•),  that  is  $  ^(<£(x))  =  x. 

The  rationale  behind  this  assumption  was  that  firstly,  was 

known  to  behave  like  2m  (see  Stephens  (1975)),  secondly,  $  ^{i/ (n+1)} 
approximates  the  i^  element  of  m  and  thirdly,  as  Vq  is  a  doubly 
stochastic  matrix  (the  sum  along  any  row  or  column  is  1)  we  may  write 


2 


W  =  r2(X,V05)  . 

De  Wet  and  Venter  (1972)  showed  that  the  asymptotic  distribution  of 

* 

W  has  the  form 

*1/2  ^ 

(2)  2n(l-W  1/^)  -  an  — >  ; 

00  2 

where  4  =  (Y^-l)/i,  {Y^,i>  1}  is  a  sequence  of  i.i.d.  N(0,1) 

variates , 

(3)  a  =  (n+1)'1  {  l  j(l-j)(«-{*_1(j)}r2  -~ 

i=l 

j  =  i/(n+l)  and  $(•)  is  the  N(0,1)  density  function. 

Beyond  the  De  Wet  and  Venter  result  the  first  step  towards  the 
asymptotic  distribution  for  W  was  to  show  that  the  Shapiro-Francia 

4* 

(1972)  statistic  W  given  by 

V*  =  r2(X,m), 

behaves  in  the  same  way  as  W  .  This  was  done  independently  and  via 

different  routes  by  Verrill  and  Johnson  (1983)  and  by  the  authors  in 

Fotopoulos,  Leslie,  and  Stephens  (1984),  henceforth  called  FLS,  where 

+  * 

expression  (2)  was  established  with  W  in  place  of  W  .  In  fact  we 
show  in  FLS  the  equivalent  result  that 

*1/2  fl/2 

(4)  n(W  -W  )  -*■  0  in  probability. 
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Our  task  in  the  present  paper  is  to  show  that 


1/2  +1/2 

(5)  n(W  -W  )  -*•  0  in  probability  . 

We  note  that  Verrill  and  Johnson  (1983)  contains  a  result  (Theorem  3) 
which  should  eventually  cover  the  asymptotic  distribution  of  W.  However 
certain  properties  of  V^m  need  to  be  established  before  it  can  be 
applied.  Inequality  (1)  does  not  appear  to  be  enough. 

2.  Asymptotic  Properties  of  W  and  a^. 

The  following  theorem  presents  one  version  of  the  asymptotic  distri- 

1/2 

bution  for  W  -  in  fact  the  asymptotic  distribution  for  W  -  whilst 
the  corollary  offers  the  complementary  form  in  terms  of  W. 

Theorem.  Under  the  hypothesis  that  the  observed  sample  is  from  a  normal 
population  the  asymptotic  distribution  of  the  Shapiro-Wilk  W  takes  the 
form: 

2n(l-W1/2)  -  2En(l-W1/2)  +  ; 

00  2 

where  C  =  (Y^-l)/i,  and  {Y^,  i  >  3}  is  a  sequence  of  i.i.d. 

N(0,1)  variables. 

1/2 

From  the  lemma  below  and  from  the  theorem  we  have  /n(l-W  )  -*•  0 
in  probability,  which  leads  to 

2n(l-W''''/2)  -  n(l-W)  =  (/n(l-W^2)  )2  -*■  0  in  probability. 

Again  applying  the  lemma  below  we  obtain, 


Corollary.  An  equivalent  form  for  the  asymptotic  distribution  of  W  is: 

D 

n(W-EW)  -*  -C  . 

It  is  not  obvious  from  their  definition  just  how  the  constants  a  will 
behave  as  n  gets  large.  The  following  lemma  should  shed  some  light  on 
this  matter. 

Lemma.  The  constants  a^  defined  in  (3)  have  the  following  properties: 

(i)  a  -  2nE{l-r (X,b) }  -+  0,  where  b  can  be  any  of  m, 
n  -  —  ■  — 

l/2V~1m  or  H, 

(ii)  a  -nE(l-W)  -  0, 
n 

(iii)  [a  -n(l-n  1m'm)  +  3/2]  <  C(logn)  \ 
n  ~  —  — 

and 

(iv)  C^loglog  (n)  <  an<  C2loglog(n),  0  *  Ci  <  C2  *  . 

Note  that  (iii)  implies  that 

m'm  =  n  -  a  -  3/2  +  0(1)  . 

-  n 

As  far  as  we  are  aware  this  property  of  m'm  has  not  appeared  elsewhere; 
the  behavior  of  m'm  is  of  interest  in  other  contexts  and  has  been  the 
subject  of  a  number  of  papers  (see  for  example,  Balakrishnan  (1984),  Ruben 
(1956)  and  Saw  and  Chow  (1966)). 
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It  should  be  pointed  out  that  the  cove rgence  for  (i)  and  (ii)  in  the 


lemma  is  extremely  slow;  for  example  a^  -  2En(l-r (X,m) )  c  -0.1  for  40  <  n 
<  400.  It  is  therefore  unclear  as  to  which  set  of  norming  constants  it  is 
best  to  use. 

When  Sarkadi  (1975)  established  the  consistency  of  the  Shapiro-Francia 
test,  it  seemed  likely  that  the  Shapiro-Wilk  test  would  share  that  property. 
That  it  is  indeed  consistent  will  follow  from  a  straightforward  application 
of  a  result  in  Sarkadi  (1981). 

3 .  Proofs . 

Notation.  We  give  some  notation  which  will  be  used  throughout  the  rest 

of  the  paper.  With  or  without  subscripts,  C  is  a  generic  constant  which 

-12  2 

is  independent  of  i  and  n.  Set  g  =  1/2V.  m,  nG  =  g'g,  nM  =  m'm, 

V  —  n  n - 

N  is  the  integer  part  of  l/2(n+l),  S2  =  E^1  (X.-X)2/n, 

nil 

^(v)  =  <*  ''■{exp(-v) } ,  si  =  I*  v  \  =  ih ( s ^ ) ,  Wi  =  -log($(X^) )-s^.  Note 


W.+s.  is 
i  i 

the 

,th  , 

l  largest 

order  statistic  in  a 

random  sample  from  an 

exponential  distribution; 

EW.  =  0,  EW?  =  d.  . 

i  i  in 

,  where  d.  =  E?  v  2, 
in  l 

EW3  =  2E? 
1  1 

-3 

and  ! EWr 

1  X  1 

<  Ci~2  for  r  >  3. 

Denote  the  element 

of  £,  m 

and 

H  by  respectively  g^,  itk  and 

H  (m  and  H  are  given 

in  section  1).  Further,  as  r  is  scale  and  location  invariant  we  assume 
without  loss  of  generality  that  our  sample  is  from  a  N(0,1)  population. 

Proof  of  Consistency.  The  consistency  of  W  follows  directly  from 
Theorem  1  of  Sarkadi  (1981) .  There  is  a  small  difficulty  in  that  whilst 
it  appears  to  be  the  case  that  V^^rn  is  a  vector  whose  elements,  as  you 
move  down  the  vector,  are  monotonic  increasing,  we  are  unbale  to  prove  it. 


1/2 

This  means  we  cannot  establish  that  W  is  always  positive.  Sarkadi 

+l/2 

exploits  the  fact  that  W  is  always  positive  to  argue  that  tests  based 

on  W  ^2  are  equivalent  to  those  based  on  W  .  We  need  to  argue  likewise 

1/2  tl/2 

for  W  (note:  we  distinguish  between  W  ,  W  etc.  and  the  square  roots 

+  1/2  2 

of  W,  W  etc.;  it  is  true  that  W  =  (W  )  but  in  view  of  what  has  just 

1/2 

been  said,  we  are  unable  to  say  whether  W  is  the  positige  square  root 
of  W) .  We  overcome  this  difficulty  by  showing  below  that 


(6) 


1/2  -1/2 

W  >  -C(logn)  ,  C  independent  of  n. 


From  the  theorem  and  the  lemma,  the  100ot%  critical  region  for  the  test 

For  the  test  based  on  W 


based  on  W1/2  is:  W1/2  <  l-l/2(c(a)+a  )n_1 

n 


it  is  W  <  l-(c(a)+a  )n 

n 


-1 


By  (6)  the  two  critical  regions  are  asymptoti- 

,1/2 


cally  equivalent.  We  need  only  show  therefore  that  W  is  consistent. 

We  extablish  (6)  by  setting  1  to  be  an  n *  1  vector  of  l’s  and  writing 


1/2  — 

U  '  =  {(X-Xl  )’(g-m)  +  X'm}/ (nS  G  )  . 

—  — n  — -  nn 


As  X'm  >  0  provided  only  that  the  components  of  X  are  increasing  (see 
Sarkadi  (1975),  Lemma  2),  and  from  (1),  maxjg^-m^l  <  C//(logn),  we  have, 
with  the  help  of  (21)  below. 


W1/2  >  -Cz"|x.xi/{nSnGnv/(logn)}  >  -C//(logn)  . 


We  turn  now  to  Theorem  1  of  Sarkadi  (1981).  Applied  to  our  context,  it 

1/2 

states  that  W  will  determine  a  consistent  test  of  H^:  that  the  random 
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sample  is  normal,  versus  H^:  that  the  observations  are  not  normal  (Sarkadi 
also  allows  the  observations  under  to  be  m-dependent  with  common 

non-normal  marginal)  providing 


(7) 


g.G'1 
l  n 


I(i-1  <  nu  <  i)4>  ■''(u)du  =  1  +  0(1) 


where  1(A)  is  the  indicator  function  for  A.  Note  that  Sarkadi's 

theorem  is  framed  in  terms  of  a  statistic  T  which  here  takes  the  form 

n 


T  =  v  {(X.-X)n  1/2S_1-c.  }2  =  2(1-W1/2)  , 
n  -  i  n  m  ’ 


where  c.  /n  =  g./G  .  To  establish  (7)  we  require  results  contained  in 
min 

the  proof  of  both  our  lemma  and  theorem,  therefore  we  will  leave  the 
derivation  of  (7)  till  the  end  of  the  article. 

Proof  of  Lemma.  We  start  by  showing  (iii);  observe  that 


n(l-M  )  =  2 v /  Var(X.)}  -  (2N-n)Var (XM) 
n  “  l  N 


We  can  write 


Var(Xi) 


E4'(si+Wi)  -  EiHs.+Wjy 


Expanding  $  in  W^  up  to  third  order  terms,  using  the  properties  of 
W^  given  in  the  section  on  notation  and  together  with  results  in  Leslie 
(1984)  (in  particular,  Lemma  6  and  the  properties  of  given  in  section  3) 
we  can  show  that 
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K'ar (X^)  -  {;■ '  (s^)  }2din  !  <  C{i(log(n/i))  }  2  , 

where  ^'(s  )  =  (expC-s^  }/ i  (i  1(exp(-si)))  and  din  =  I?  v  2.  This  yields 
N  N 

(8)  1 1  Var  (X . )  -  :  {«'(s.)rd.  |  <  C(logn)"z  . 

1  1  1 

Using  the  Euler-Maclaurin  summation  formula  (Knopp  (1951),  p.  534) 

(9)  0's  -log((n+l)/i)  -  l/2(r1-(n+l)_1)  <  {i~2-(n+l)"2}/12 

and 

(10)  0 <  d  -  i_1{l-(i/(n+l))}  -  l/2{i_2-(n+l)~2)  <  l/(6i3)  . 

In  FLS  we  show  that  ! v  * (v) j  and  |v"(v)l  are  monotonic  decreasing  in  v; 
also  in  Lemmas  1  and  4  in  Leslie  (1984)  it  is  shown  that 

(11)  '  {log((n+l)/i)  }  \  <  C { log (n/ i) }  ^ 
and 

(12)  | log(  (n+1)  /i) }  |  <  C{log(n/i)}  ^ 2  . 

With  (9),  (10)  and  (11)  we  have 

(13)  ! {v ’ (s . ) } 2 (d  — i_ 1 { 1— ( i / (n+1)) )) ]  <  Ci~2/log(n/i) 

i  n 


9 


r 


(14)  )  i-v'  (s.)}2  -  (4 '  {log(  (n+1)  /  i ) } ) 2  {  <  C |y ' (a . )y"  (a . ) j /i  , 

where  log{(n+l)/i}  <  <  s^. 

Expressions  (11)-(14)  taken  together  imply  that 

N  N 

(15)  |E  {*'(s.)}2d.  -  (n+l)~2  l  (?{<r1(i/(n+l))})"2i{l-i(n+ir1}| 

1  ir  T  1  in  ^  1 

<  C(logn)  1  . 

From  the  definition  of  a  and  with  (8)  and  (15)  we  obtain  (iii) . 

n 

Next  we  establish  (iv) .  A  well  known  inequality  is  useful  here  (see 
Renyi  (1970)  p.  164;  for  x  <  0 

(16)  <f(x)(l-x  2)/jxj  <  0(x)  <  (f(x)/|xj  . 


From  this  we  obtain,  for  1  <  i  <  N,  and  with  x  =  H., 

=  1 


(17) 


1-H . 2  <  i;H. I / { (n+l)b (H . ) }  <  1 


1  =  '  1 


In  view  of  the  symmetry  in  the  summands  in  a^,  we  need  consider  only 
1  <  i  <  K.  We  use  (17),  over  the  range  1  <  i  <  [-|-N]  and  for 
[yN]  <  i  <  N  we  use 


(18) 


C.  <  <)> (H. )  (i  / (n+1) )  { 1—  ( i /  (n+1) ) }  <  C0  , 

i.  1  A. 


where  C^,  do  not  depend  on  i  or  n.  Based  on  (16)  we  show  in  Lemma  3 
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of  FLS  that  for  any  cQ(0  <  < 1)  there  is  a  y(cq)  such  that  when 

0  <  u  <  Y (Cq)  <  -  , 

(19)  -{-log(27TU2)}1/2  <  $'1(u)  <  -{-c0log(2TTU2)}1/2  . 


I 


This  yields  for  1  <  i  <  N, 

(20)  C3(log(n/i)}1/2  <  JH± 1  <  C4{log(n/i))1/2  . 

Applying  (17),  (18)  and  (20)  we  find 

[ifil 

C5  l  {ilog (n/i) }  +  C6  <  an  +  3/2  <  C?  £  {ilog(n/i)>  +  Cg 

which,  after  approximating  the  sum  by  an  integral,  establishes  (iv) . 

To  complete  the  lemma  we  prove  (i)  and  (ii) .  First  however,  we  need 
two  results  which  will  be  used  here  and  in  the  proof  of  the  theorem: 

(21)  Gn  -  i  as  n  -*■  00 ,  and 

(22)  0  <  II  mil  //  g//  -  m'£  <  M  G_1//  £-m//  2  . 

—  —  =  n  n  ~ 

It  is  well  known  that  M  -*■  1  as  n  -*•  «  (see  Hoeffding  (1953)).  On 

22  -1  2-1 
writing  =  +  2m'(£-m)n  +  //  g-m//  n  ,  from  (1)  and  Schwarz 

inequality  we  obtain  (21) .  We  demonstrate  (22)  by  exploiting  an  idea 

in  Sarkadi  (1972).  First  note  that  m'£  >  0,  for  m'£  =  m’V  ^in  and 

being  a  covariance  matrix,  is  positive  definite.  Set  9  to  be  the  angle 


1 


between  m  and  £  then  cos  6  >  0  and  0  <  9  <  — rr.  Consider  the 

triangle  formed  by  vectors  m,  £  and  a  =  m-j>,  respectively  lines  AB, 

AC  and  CB.  Let  CD  be  the  perpendicular  from  C  to  AB.  Then 

//  a//  2  >  (CD)2  =  //  g//  2sin20  =  //  g//  2(l+cos6)  (l-cos9)  >  // g//  2  (l-cos6)  >  0. 

As  cos9  =  m’g{//m//  //  g//  },  (22)  follows. 

Returning  to  the  proof  of  (i)  and  (ii)  of  the  lemma  we  show  first 

that 

(23)  |nE(l-r(X,b))  -  n(l-Mn)  +  3/4 |  <  C(logn)-1  . 

2 

As  r  is  scale  invariant  and  as  is  sufficient  for  the  scale  para¬ 

meter  a  we  can  use  Theorem  7,  p.  243  of  Hogg  and  Craig  (1970)  to  yield 

nE  r  (X,b)  =  m'b/{ESn//  b//  n'1/2}  . 

2  2 

With  nS^  distributed  as  x  on  n-1  degrees  of  freedom  it  is  elementary 
to  show  that 

ESn  =  (2/n)1/2r(n/2)/F((n-l)/2)  . 

-1  -2 

By  Stirling’s  formula  this  reduces  to  l-(3/4)n  +  0(n  ).  As 

-1/2 

n  //  b //  -*■  1  (the  case  Id  =  H  is  shown  in  Lemma  2  of  De  Wet  and 

Venter  (1972)),  and  using  (1),  (22)  and  an  analogue  of  (22)  with  £  and 

G  replaced  by  H  and  H  =  /{(H'H)/n}  (this  analogue  holds  because 
nr  —  n  - 

m^H^  >  0  for  all  i,  nu  and  always  having  the  same  sign)  we  have, 


12 


En{  1-r  (X,b) }  =  n(ESn)_1{l-(3/4)n_1  -  n“1/2m'b// b// _1}  +  0(n_1) 

=  n(ES  )_1 (1-M  )  -  (3/4)  +  O(logn)'1  , 
n  n 

*  n(l-Mn)  -  (3/4)  +  O(logn)-1  , 

the  latter  expression  resulting  from  the  fact  that  n(l-M^)  =  O(loglogn) 

(using  (iii)  and  (iv)  of  the  lemma  and  recall  that  M  -*•  1) .  This 

n 

establishes  (23).  Analogous  to  (23)  for  b  =  £  we  have 

(24)  |nE{l-r2(X,£)}  -  n(l-M2)  +  3/2 |  <  C(logn)-1  . 

2 

To  show  this  we  note  that  as  nES  =  n-1  we  can  write 

n 

nEr2(X,£)  =  E(X'^)2/{(n-l)G2} 

with 

e(x'£)2  =  £'voi  +  te's)2  =  ^b's.  +  (^’e.)2 

=  4nM  G  +  (n  G  M  )2  +  0(n/logn) 

2  n  n  n  n 

using  (1)  and  (22).  Again  using  the  property  that  n(l-M^)  =  O(loglogn) , 
we  obtain  (24) .  As 

(25)  2n(l-M  )  -  n(l-M2)  =  {/n(l-M  )}2  =  0{ (loglogn) 2/n} 

n  n  n 

it  is  clear  from  (23),  (24)  and  (iii)  of  the  lemma  that  (i)  and  (ii)  hold. 
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2 

Undoubtedly  it  is  true  that  a  -E^d-r  (X>b.) }  0,  for  b_  =  m 

and  H.  However  this  entails  showing  that  //  VqH-2H//  -*•  0  and 
//  V^m-dm//  -*•  0  both  of  which  will  follow  once  is  replaced  by  the 

approximation  V  given  in  Leslie  (1984):  corollary  1  in  Leslie  (1984) 
permitting  this.  These  two  results  will  involve  a  quantity  of  tedious 
analysis  and  it  seems  unnecessary  to  set  it  down  here. 

Proof  of  Theorem.  The  theorem  follows  from  (2),  (4)  and  (5)  together 
with  the  lemma;  therefore  to  prove  the  theorem  it  remains  to  establish 
(5).  Now 


nS  (W1/2-W+1/2)  =  l  X. (g.G'^-m.M-1) 
n  “  i  °i  n  i  n 


n 


,-l  ..-1> 


=  I  (Xi-mi)  (gi-mi)Gn1  +  I  (Xi-mi)m1(G'J'-M”'1')  +  (m’£ -  //  m//  //  g//  )G^ 


-1  ..-1, 


-1 


As  S  -*•  1  a.s.  and  with  (21)  and  (22),  expression  (5)  will  follow  from 
n 

Markov's  inequality  once  we  demonstrate  that 


(26) 

Ejz(X.-mi)(g.-mi)  |  -* 

(27) 

E  j I (X . -m . ) (G_1-M_1) 1 
ii  n  n 

Result  (26)  follows  from  Schwarz  inequality: 


E  j  Z  (X.-m. )  (g  ,-m.  )  j  <  {n(l-M2) //  £-m//  . 
liii  -  n  — 

With  (1)  and  with  (iii)  and  (iv)  of  the  lemma  we  have  (26). 
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To  deal  with  (27)  we  note  that  in  Lemma  11  of  FLS  we  show 


(28) 


E ! X i  — d-1  i  I  <  C//{  ilog  (n/ i)  } 


and  in  Theorem  1  in  FLS  we  show  that 


(29)  <  Ci_1(log(n/i) }  3j/2  ; 

both  of  these  bounds  hold  provided  1  <  i  <  N.  As  |  //in//-//  g//|  <  //m-g//  , 


(30) 


1  n  n 


<  //m-g//  /  (M  G  /n) 
=  —  n  n 


C(nlogn) 


-1/2 


and 


(31) 


E|Z(Xi-mi)mi!  <  2  [  (E lx±-«± |  |mj[  |  |tt>i | ) 


From  (29),  (9),  (20)  and  the  monotonieity  (decreasing)  of  |i|/(v) 


(32) 


| m  |  <  C(log(n/i)}1/2,  1  <  i  <  N  , 


so  by  combining  results  (28)  to  (32)  we  find 


! I (X.-m. )m. (G  <  C(logn)_1/2 

l  l  i  n  n  = 


This  establishes  (27)  and  hence  the  Theorem. 
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Derivation  of  Expression  (7).  Denote  the  integral  in  (7)  by  J(i,n)  then 


J(i,n) 


1(l/n)},  for  i  =  1  , 

<  4>  ^((i-l)/n)n  ^  +  -jn  2(i£{i>  ^  ( (i-0)  /n) })  \  0  <  6  <  1,  l<i<n 
${$  "'■(l-n  1)}  ,  for  i  =  n  . 


Without  loss  of  generality,  assume  n  is  even.  Then 

1 

n  2 n 

I  g.J(i,n)/G  =  2n  1  £  g .  (4>  1{  (i-l)/n}  +4n_1  {${$  1((i-9)/n)}}  1) 

1  1  n  i=2  2 

+  2gn<J> { 4>  1(l/n)}  . 


By  (16),  for  1  <  i  <  -|n. 


${$_1((i-0)/n)})  >  (i-1) /n)  > 


(i-l)$  1{(i-l)/n}/n,  2<i<kn,  k<~ 
C(k) ,  kn  <  i  <  yn  . 


Thus  by  Schwarz  inequality, 

1  i  I 

2  n  1  len  2.  2  ^ 

|n  2  £  g./${$  ^  ( (i— 6 )  /n) }  j  <  n  2G  {  £  (i$  ^(i/n)}  2  +  £  C(k)  2n  2) 

2  1  n  1  kn+1 

-1/2 

which  in  turn  is  bounded  by  C{nlog(n)}  ,  in  view  of  (19).  Further,  by 

(16)  and  (19),  <K$  ^(1/n))  'v  0((logn)/n),  by  (1),  gn  ^  mn  and  with  (32) 
and  finally  (22)  we  can  argue  that 
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1 

_ -j  ^  —1  —1  ^  —1  1 

2n  I  g.$  {(i-l)/n}  ^  n  7  g.<5>  (i/(n+l)}  ^  n  m'g  'v  M  G 
2  1  1  1  —  A  n  n 

These  ensure  that  (7)  holds. 
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